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loss happening during the operation. One of the major causes is cold air 
leaking from the inside of the cold room to outside. It is due to poor 
Keywords: insulation and improper maintenance of cold storage. It is very difficult to 
identify the exact location of the leak by the naked eye. In this research work, 
the cold storage was inspected with the thermal imaging system. Thermal 
l cameras are highly sensitive to temperature and it can detect the variation of 
Insulation problem temperature ranging from 0.1°C. The measured temperatures are further 
Thermal image converted into a colour based pattern. It is known as thermogram. These 
colour-based thermal patterns are further processed for identifying energy 
transmission location. It is done by applying various image processing 
methods such as histogram equalization, diffusion error, otsu thresholding 
and morphologic function. These techniques were applied to images of cold 
storages and exact cold air transmission locations were identified. 
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1. INTRODUCTION 

Fruits and vegetables start to deteriorate or lose their quality in the following ways: change in 
weight, texture, flavour, colour, nutritional value and appearance at post-harvesting temperature and tends to 
reduce the shelf life of the product. Deterioration of the products depends on storage time and 
temperature [1]. Post-harvest cooling has the potential to reduce the respiratory activity, internal water loss 
and wilting of the products. It is also helpful to prevent degradation by enzymes, stops the production of 
a ripening agent, ethylene and the growth of microorganisms causing decay and pathogens. Moreover, post- 
harvest cooling has the tendency to maintain the quality of the products for an extended period and supply it 
whenever required. It is the important benefit for growers and distributors. Recent day cold storages are 
highly capable of maintaining organoleptic qualities of the products and reduce spoilage and increase 
the shelf life [2]. Infiltration of water moist air into cold storage unit cause many problems which includes 
increased cost of running, defrosting the refrigeration system, mist formation [3]. Moreover, the temperature 
fluctuation causes weight loss and quality deterioration of the products [4]. Prediction and measurement of 
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the cold air infiltration through doors by an analytical model has been proved in that Gosney model is better 
for predicting infiltration [5]. Cold storage consumes more energy than common buildings for maintaining 
lower temperatures. In the food industry, the refrigeration system is the only technology with the high cost 
and produce more emission, its approximately 1 % of global GXG emissions [6]. It very well may be seen 
that the energy demand varies significantly from various kinds and areas of cold storages. Meanwhile, cold 
storages are the structures with huge warm mass, the transient energy utilization of which much depends on 
the variety of cooling heap of interior mass. To full fill the requirement, it is the necessity of cold storage 
exceptionally relies on its types and area. Hence, with an enormous energy devouring per storage volume, for 
the reasons to discover retrofit exponential of existing cold storage and better handling of the energy supply, 
an examination and quantization of the influencing variables of the transient energy consumption is important 
for users [7]. The innovation of phase change materials utilized in the cold room has incredible potential in 
energy storage and cost sparing under the foundation of peak and valley price of electricity [8]. 

In recent days the thermal imaging system has been used as an inspection tool for measuring 
different temperatures. It is working with the principle of infrared radiation (IR) to absorb the temperature 
and the radian energy of the object. Its major advantages are non-contact and not alter the temperature of 
the measuring objects. Moreover, it has the potential to measure wide temperature ranges, highly sensitive, 
low power consumption, and long life. In addition to that it is easy to operate, reliable and safe [9]. In Non- 
destructive testing, TI has been used to detect the material's defects and inhomogeneity by its strange thermal 
behavior [10]. Identification of a defect in surface and under the surface can evaluate by inspecting 
the contrast of the thermal image pixels. It is done by using the undamaged area pixel contrast 
as reference [11]. In civil engineering TI system is used to inspect and ensure the quality of the building by 
detecting the surface defects in structures [12]. In this research work, cold rooms are inspected with thermal 
camera and images are processed and analyzed with image processing tool in MAT lab. 


2. MATERIAL AND METHODS 

In this experimentation, two-ton capacity cold storage was used. The cold storage has been fully 
loaded with potatoes, which was collected from local agricultural lands in Chennai. Potatoes are stacked in 
the plastic bins and then it was placed in order. As per the regulatory guidance for storing temperature of 
potatoes, the room temperature was maintained 10°C with 90% humidity range [13]. The inside temperature 
was controlled by a control panel placed outside the cold storage. The environment temperature of Cold 
storage has also been maintained with a consistent temperature. The output and input parameters of the cold 
storage were regularly monitored and noted. A detailed description of the cold storage unit is provided in 
Table 1. To inspect the cold storage, FLIR E60 model camera used. Before taking the images, the outside 
environment of the cold storage was set as constant temperature. The automatic control system has been used 
to change ambient temperature and relative humidity of the cold storage. The calculation of the transmission 
heat created by walls, floor and ceiling requires information on thickness and type of isolation material used 
in construction of cold room, construction of building, physical specifications of the cold storage volume, 
inside and outside environment temperatures. Following equations were used to calculate transmission 
heat [14, 15]. 


qc = KA (tout =f) 

Qc = transmission heat at flat surface (kCal) 

K= Total heat transmission coefficient (kCal/hm?°C) 
A= Area of heat transmission (m7) 

to= Temperature of outside or neighbour volume (°C) 
ti= Inside volume temperature (°C) 


a_i n Xd 
K ai + 1=1" jj + ao (1) 
Where: 


Oo= Coefficient of heat transmission of outside surface (kCal/hm? °C) 

ai= Coefficient of heat transmission of inside surface (kCal/hm? °C) 

xi = Material thickness (m) 

i = Thermal conductivity (kCal/hm °C) 

The proposed method helps evaluate a leak of cold air through the introduction of image processing 
techniques. The following steps were taken to evaluate the exact location of the leak. 

1. Acquisition of the image using a thermal camera. 

2. Screening the image using filters. 

3. Enhance the image. 
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4. Compare the results 








Figure 1. Flow chart the system 


The image processing techniques currently evaluated are histogram, Otsu thresholding, diffusion 
error and morphological function (1) show the flow chart of the system. The purpose of this research is to 
detect a cold air leak location in a cold storage unit. 


3. FLOWCHART AND EXPLANATION 
3.1. Image acquisition 

The quality of the thermal images depends on controllable parameters such as ambient temperature, 
airflow and lighting. It is necessary to maintain the temperature at 30°C, without air draft and diffused 
lighting. Thermal camera has been used to inspect the cold storage after 24 hrs of installation. The thermal 
camera switched to video mode and recorded the entire cold storage from the outside. It is observed that 
the thermal camera shows most of the region with yellow colour and some of the regions with a blue colour. 
Based on this variation the blue spot region is marked and images of that region were taken for further 
analysis. Thermal images were captured at where the low-temperature region is identified. Each image is 
analyzed and few images are selected for image processing. 


3.2. Histogram method 

Two sets of images are considered for the analysis-thermal cold storage image and trial image for 
training. The trial thermal pattern was generated using computational fluid dynamics software (Star CCM+). 
The trial image is generated such that, the colour blue depicts the lowest temperature and extends to 
the colour red which stands for the highest temperature [16]. The analysis was performed in the numerically 
generated image and then repeated for the thermal image of the cold storage 


Table 1. The cold storage unit 


Image Name Zeros Non-Zeros RMS Error SNR Median Stddv 

Imagel (No leak) 52 204 0 0.000 206 9.7e+002 
Image2 (Moderate leak) 0 218 65.37 -9.94 193 4.3e+002 
Image3 (Light Leak) 38 256 54.01 -7.53 185 5.5e+002 
Image4 (SevereLeak) 0 256 80.54 -11.64 190 2.5e+002 
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3.2.1. RGB to grey conversion 
RGB is converted into a greyscale image, which is the intensity image. This conversion occurs by 
eliminating hue and saturation and retaining luminance. 


3.2.2. Noise filtering 

The median filter is used in this paper for sorting out salt and pepper noise. It keeps the sharp edges 
and is a non-linear operation. For every input pixel p(x, y), the value of the pixel is found, its median is 
identified and the value is assigned to output pixel q(x, y). After acquiring the image, it is processed using 
the three types of processing techniques, namely histogram, Otsu thresholding and morphological function. 


3.3. Otsu thresholding method 
The acquired image is first converted in to a grey image. The red component is extracted from 
the image and filtered using a median filter. 


3.3.1. Median filter 

The median filter is used in this method for preserving the edges and to remove the salt and pepper 
noise. The process of the median filter is the replacing of each value with the median value of 
the neighbouring pixel. It is determined by first arranging all the pixel esteems from the window into 
numerical request, and after that supplanting the pixel being considered with the centre pixel esteem. 
The focal pixel is supplanted by the middle of pixels in the window. 


3.3.2. Histogram equalization 


The image is processed using histogram equalization which helps adjust the image intensities to 
enhance contrast. 


Table 2. Histogram analysis 


No Image Histogram value Observation 

1 Image 1 No peaks in the ranges less than 150 Side wall 

2 ĮImage2 Medium pixel values in <150 region Moderate leak 
3 Image 3 less pixel value in <150 region Light leak 


4 Image 4 More pixel values (peaks) in <150 region severe leak 


3.3.3. Otsu thresholding 

In image binarisation, Otsu threshold method helps to reduce the grey level image to a binary image. 
In this method, all possible threshold values are iterated to calculate the spread of pixel value on both sides of 
the threshold value. The purpose is to find out the threshold value for classifying the is either leaking or non- 
leaking. According to OTSU method, minimizing ‘within-class variance’ is equal to maximizing “between 
class variance' of the segmented image. However, maximizing ‘between class variance’ is less expensive than 
minimizing ‘within-class variance’. 


3.4. Morphological process 
3.4.1. Gabor filters 

The RGB image is converted to a grey image. The captured image is filtered using a Gabor filter. 
The Gabor filter converts the image to discrete value for getting the constant region of the image. Compared 
to other filters, the Gabor filter is primarily used for a low range of values. The filtered image is iterated for 
five iterations. This helps extract particular affected area for a clearer view. 


3.4.2. Enhancement techniques 

In this technique, there are two types - dilation and erosion. Dilation is the process which increases 
the white content of the image by adding pixels to the boundary of the object. Erosion is used to increase 
the black content of the image by removing the pixels from the boundary of the object. The small particles in 
the dilated image are eroded in the erosion process. 


3.4.3. Morphological function 

A morphological method identifies the shape and features of a cold air leaking area in a given 
image. The organizing component is set in every conceivable situation in the picture and it's contrasted and 
the relating neighbourhood of pixels. According to the output of the test, it identifies the image that fits or 
hits the neighbourhood pixels. Erosion strips one layer of both inner and outer boundaries. It avoids small 
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details and enlarges holes and gaps between different regions. The boundaries are found by subtracting 
the eroded image from the original image. x=f—(f(@)s), where f is an image of the regions, s is a 3x3 
structuring element, and x is an image of the region boundaries [16]. Dilation adds one layer on both 
the inner and outer boundaries. It will fill in the gaps between the two regions. 


4. RESULT AND DISCUSSION 
The analysis is done with four images taken at different location of cold storage. Each image is 
taken separately and the results are observed using four methods and compared. 


4.1. Histogram method 

According to Figure 2, it was observed that the pixel intensity of image 1 (a) lies between above 150 
to 255 range. It is due to a high number of bright sides by high temperatures. It is noted as no leak of cold air. 
At the same time pixel intensity of image 2 (b) has shown a light peak in the histogram at the range less than 
150. It is due to darker regions of images caused by low temperatures. Similarly, the histogram of image 4 (d) 
shows more peaks in the lower range and image 3 (c) has a moderate peak in lower ranges. It is observed that 
the histograms show the cold leaking regions clearly by the distribution of pixel intensities. 








Figure 2. Histogram of the images 
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4.2. Otsu thresholding method 

According to Figure 3, the same images which has been used in the histogram method are processed 
with otsu thresholding. The result shows the image 4 (d) having more leaking regions than others. 
Image 1 (a) shows no leaking, similarly the images 2 (b) and image 3 (c) showing less and moderate 
leaking respectively. 


Moderate leak 





(c) 


Blue s 





KINAI- 


(d) 


Figure 3. Thresholding method of images 


4.3. Error diffusion dither 

Same images were also used in this thresholding. After histogram equalization of all the images, 
error diffusion dither thresholding applied to all the images and the results are comparable. In this method, 
the exact location of the leak was identified. Image 4 Figure 4. shows more leaking region than other images. 
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Figure 4. Error diffusion dither 


4.4. Morphological method 

According to Figure 5, the thermal image (a) is converted into a greyscale image, and the intensity 
of the image is noted (b). The grey image has undergone five iterations and all the stages are shown in (c). 
After iteration, the Fig.ure proceeds with the morphological process. For extracting the leak area, the iterated 
image analyses with the morphological process (dilation and erosion) (d). The leak location of the image is 
extracted and it is incorporated with the original image. This helps to identify the leak location without an 
expert's advice (e). 


ia 
128 


Figure 5. Morphological functions 








5. CONCLUSION 

Thermal images can be used to identify the cold air leak location of cold storage by analyzing 
the images, which helps to predict early, the cold air leakages and prevent further energy loss. Four methods 
have been used for the thermal image analysis. It is observed that all four methods are useful. In that, 
the histogram method helps identifying a leak location. However, the other three methods help in extracting 
the region of interest by enhancing the leak location. In future, it is recommended that the analysis has to be 
done by taking images from inside the cold storage and finding the heat transmission into cold storage from 
outside. Moreover, image classification algorithm should be made for detailed analysis and improve 
the accuracy of results. 
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